Mon. Not. R. Astron. Soc. OOP, [Tp2l (2008) Printed 18 March 2009 (MN I^TeX style file v2.2) 



Fragmentation in Molecular Clouds and its connection to 
the IMF 



Rowan J. Smith^ Paul C. Clark^ and Ian A. Bonnell^ 

^ SUPA, School of Physics & Astronomy, University of St Andrews, North Haugh, St Andrews, Fife, KY16 9SS, UK 

^ Zentrum fiir Astronomic der Universit at Heidelberg, Institut fiir Theoretische Astrophysik, Albert- Ueberle-Str. 2, Heidelberg, Germany 



18 March 2009 



ABSTRACT 

We present an analysis of star-forming gas cores in an SPH simulation of a Giant 
Molecular Cloud. We identify cores using their deep potential wells. This yields a 
smoother distribution with clearer boundaries than density. Additionally, this gives an 
indication of future collapse, as bound potential cores (p-cores) represent the earliest 
stages of fragmentation in molecular clouds. We find that the mass function of the 
p-cores resembles the stellar IMF and the observed clump mass function, although 
p-core masses (~ 0.7 M©) are smaller than typical density clumps. The bound p-cores 
are generally subsonic, have internal substructure, and are only quasi-spherical. We see 
no evidence of massive bound cores supported by turbulence. We trace the evolution 
of the p-cores forward in time, and investigate the connection between the original 
p-core mass and the stellar mass that formed from it. We find that there is a poor 
correlation, with considerable scatter suggesting accretion onto the core is dependent 
on more factors than just the initial core mass. During the accretion process the p- 
cores accrete from beyond the region first bound, highlighting the importance of the 
core environment to its subsequent evolution. 
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1 INTRODUCTION 

The origin of the stellar initial mass function (IMF) has been 
a major question in star formation since it was first measured 
by [Salpeter (1955). The similarities between the IMF and 
the mass function of clumpy structure in molecular clouds 
(MC's), has led many to propose a link between the two (e.g. 
[Motte et al.|1 998'). In this scenario, stars are formed directly 
from the dense cores of gas observed in molecular clouds 



(e.g. [Johnstone et al.|2000l|Testi Sargent] 1998) [Johnstone 



et aL||2006| [Nutter &; Ward-Thompson[|2007| [Enoch et al , 
2008| . A direct link between the core mass and the resulting 
stellar systems mass is often assumed. For example [Alves| 
et al. ( 2007 ) propose an efficiency of one third between core 



mass and stellar system mass. 

This fits into the the long-standing idea that fragmenta- 
tion of molecular cloud structure forms the IMF. This could 
be simply due to gravity ( Larson|1973 Elmegreen & Math- 
leu [19831 , causing successive fragmentation of a larger body 



as it collapses. Alternatively, it has been shown that the 
supersonic turbulence observed in MC's produces a hierar- 
chical density structure, the dense peaks of which have a 
mass distribution similar to the IMF ( Henriks'eii|[1986 Lar- 
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son IW; 'Elmegreen"1997"; •Klessen"200l] [Padoan Nord^ 
fund 2002,i Hennebelle Chabrier 2QQ8| . The Jeans mass 
^ Jeans||1902 ) at the point of fragmentation has been shown 
to be only weakly dependent on temperature, density, met- 
alicity and radiation field in the environments in which stars 
form ( [Larson||20Q5| [Elmegreen et al.||2QQ8| , which leads to 
the characteristic mass of a Jeans unstable fragment being 
similar in all molecular clouds. 

However, it is unclear how to get to the IMF from the 
clump mass function. Many assume a direct 1 — 1 link be- 
tween the masses (Motte et al. 1998[ [Padoan & Nordlund 



20Q2[ [Alves et al.|[2007[), while ot hers include the effects of 



multiplicity ( Goodwin et al.|2008| ). However, there are many 
complicating factors in this story of collapse, such as feed- 
back from winds and outffows (Shu et al. 1988 Silk 
Myers 2008), supporting magnetic ffelds (Heitsch et al 



Tilley & Pudritz 2007 ) and competitive accretion ( Zinnecker 
1982( [Bonnell & Bate[[2006[). AH of these processes are m- 
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2001 



volved in the collapse of a fragment to a star and all could 
vary locally. In fact, Swift & Williams (2008) have shown 



that when a core mass function is evolved into a stellar IMF, 
a Salpeter like distribution was found regardless of whether 
the core-to-star efficiency was constant, variable or included 
multiplicity. Moreover, under the competitive accretion the- 
ory of star formation there is no need for a direct correlation 
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between core masses and stellar masses at all, as the cores 
can be thought of as 'seeds' from which accretion will build 
up the future IMF (IClark & Bonnen""2005"). 



Resolving the issue of core evolution is further compli- 
cated as their lifetimes are not well known. In the classi- 
cal Shu (1977) picture a quasi-static core supported by its 
magnetic field will slowly collapse to form a star. However, 
in the more dynamical view, usually proposed in gravotur- 
bulent fragmentation, cores can collapse quickly when they 
become unstable, making the process hard to observe. Fur- 
ther, the variation of core free-fall time with density means 
that the clump mass function observed at a snap shot in 



time might actually evolve into a steeper IMF ( Clark et al.| 



2007 ). 

Massive star formation poses another question for the 
evolution of cores. Any core which is large enough to form 
a massive star directly is likely to fragment without some 
additional heating mechanism. McKee & Tan (2003) the- 



orize that massive cores could be supported by their in- 
ternal turbulent energy, although simulations have shown 
that some fragmentation is unavoidable (Krumholz et al. 
Dobbs et al.,,20Q5t . Alternatively competitive accre- 
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tion predicts that there will be no massive cores that do not 
fragment into smaller structures. It would be one of these 
smaller cores which would preferentially accrete from its en- 
vironment to become a massive star (Bonn ell et al.|[2004 ). 
In this case, the link between core mass and stellar mass is 
destroyed by accretion. 

Observationally it is hard to be sure which structures 
will gravitationally decouple from their environment and are 
hence 'pre-stellar' in nature (i.e. will form stars in their fu- 
ture). Observations of molecular tracers such as CO produce 
a core mass function resembling the IMF ( e.g. |Ikeda et al 



2007|. However, 
est gas. AsjLada ( 
takes place above 



these tracers are insensitive to the dens- 
1992) has shown, star formation generally 
a density of n > lO'^ cm~^. Observations 



of CO cores are typically more massive and on larger scales ( 
e.g. Tachihara et al.|2002 ) than mm-continuum cores which 



trace denser gas, ( e.g. Motte et al.||1998 Johnstone et al 



[2000 ) . Moreover, it is often difficult to define a 'core' from 
the data, without invoking a somewhat arbitrary boundary 
(e.g. [Padoan et aL 2006| Pineda et al. 2008 Schnee et al. 



2008}. In particular. Smith, Clark Bonnell (2008) demon- 



strated that core properties are extremely sensitive to the 
core boundaries, which in turn depend on the resolution, 
density range and dimensionality of the dataset. High reso- 
lution observations of dense gas in nearby molecular clouds 
are the most likely to find pre-stellar objects, for example 
those identified from a synthesis and re-analysis of the liter- 



ature in Ophiuchus by Simpson et al. (2008). However, these 



observations still suffer from a lack of completeness, and of 
course all observations are necessarily at a snap-shot in time 
and identify cores in a variety of evolutionary states with no 
guarantee that they will form stars. 

In this work we seek to determine the relation between 
the properties of the cores formed through the fragmentation 
of molecular clouds and the 'stars' which form from them. In 
our SPH simulation of a molecular cloud we identify bound 
pre-stellar cores with well defined boundaries in a complete 
data-set without time effects. Our cores are gravitationally 
bound, and as they cannot be temporarily confined by pres- 



Table 1. Given below are the initial conditions of the simula- 
tion analysed in this paper. The mass resolution is the minimum 
mass gravitational forces can be resolved for and is calculated via 

Mres - imMtotal/Npart. 



Size 

Mass 

Particles 

Mass resolution 

Dynamical time 



3 X 3 X 10 pc 
IC^ Mo 
5.5 X 10^ 
0.18 M© 

4.74 X 10^ yrs 



evolution from when first bound through their early evolu- 
tion towards stars. Section [2] outlines the simulation initial 
conditions, and Section |3] discusses how we define our cores, 
here called p-cores, using their gravitational potential to en- 
sure they are bound with respect to their environment. In 
Section |4] we calculate the p-core properties, compare them 
with observations and find the core mass function. In Section 
|5]we link the p-core masses when they are first bound with 
the mass which will be accreted by the stars (sink particles). 
Finally in Section [6] we discuss our results. 



2 THE SIMULATION 

We use a three dimensional smooth particle hydrodynamics 
(SPH) code (Monaghan 1992), to simulate star formation in 
a collapsing gas cloud. The SPH code has variable smoothing 
lengths, artificial viscosity ( Gingold k, Monaghan|1983 ) and 
pressure gradients. Self- gravity is calculated via a binary 

(Bate et aI][T995t 



tree (Benz 1990) and sink particles (liate et ai.||lljU 5|) are 
used to model the star formation and prevent the simulation 
exceeding its mass resolution. 

The initial conditions consisted of a cylinder containing 
10^ M© concentrated at one end so the top was over-bound 
and the bottom under-bound. A turbulent velocity grid, con- 
sistent with a Larson velocity dispersion of a oc r°"^, was 



generated according to Dubinski et al. (1995) and Myers 



& Gammie ( 1999 ) and interpolated onto the SPH particles. 
The magnitude of the turbulent velocities is chosen such 
that globally the cloud is initially supported by turbulence, 
which equates to a r.m.s. velocity of 4.7 kms~^. The turbu- 
lence is not driven but additional kinetic energy is released 
by the gravitational collapse of the MC, when the simula- 
tion was terminated the r.m.s. velocity had decayed to 3.65 
kms""*^. 

Table [l] shows the properties of the simulation. The 
cloud is modelled with 5.5 x 10^ SPH particles, which gives 
a mass resolution of 0.18 M© ( [Bate Burkert| |1997). The 
simulation was run on the SUPA Altix computer at The 
University of St Andrews. 

A barotropic equation of state is used for basic heating 
and cooling to ensure that the Jeans mass at the point of 
fragmentation matches the characteristic stellar mass. 



where 



sure they will collapse (c.f. Dib et al. 2007 ). We trace the core 
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Figure 1. The simulated Giant Molecular Cloud. The colours 
represent column densities in the range O.OOlgcm"^ (red) to 10 
gcm~^ (white). Sink particles are shown as white dots and the 
cloud is viewed along its long axis. 



and pi = 5.5 x 10~^^ gcm"^,p2 = 5.5 x 10~^^ gcm"^ and 
P3 = 2 X 10"^^ gcm"^ 

This equation of state mimics the effects of line cooling 



(Larson 2005 Jappsen et al. 2005) and then dust cooling 



when the dust is coupled to the gas ( Masunaga & Inutsuka 



2000 ) . When the gas becomes optically thick to IR radiation 



the gas will heat again. We invoke this heating at somewhat 
earlier stage in the collapse than typical (p2 — 5.5 x 10~^^ 
gcm"*^) to ensure the Jeans mass of a fragment is always 
resolved. Finally, we return to an isothermal equation of 
state at y93 = 2 X 10~^^ gcm~^, to allow sink particles to 
form, as we require the fragments from which they form to 
be both bound and collapsing. Therefore in our simulation, 
we will not create any sink particles until at least a density 
of p3 = 2 X 10~^^ gcm~^ has been reached. Sink particles 
will accrete bound material within a radius of 200 AU, and 
have their mutual gravitational interactions smoothed to 40 
AU. 

The simulation evolves with self gravity until just after 
a dynamical time it is as shown in Figure^ The simulation 
is terminated at 1.4 idyn^ at which point there are 949 sink 
particles which have a combined mass of 1249 M©. This 
equates to a star formation efficiency of 12.5%. 



3 CLUMP FINDING USING POTENTIAL 



As previously shown in Smith et al. ( 2008 ), clump mass func- 
tions always retain a universal profile, but the sizes, masses 
and positions of cores vary with the method of extraction. 
This behaviour is also seen in observations. IMotte et al.l 



(1998) and Johnstone et al. (2000) in their observations of 



p Ophiuchus find similar clump mass functions but the po- 
sitions of their clumps and their masses do not correspond 
with each other. As the density distribution of molecular 
clouds is hierarchical with structure on all scales it is un- 
clear which scale is the most relevant for star formation. 

In this study, we take a different approach, since we are 
interested in establishing the connection between the mass 
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Figure 2. A cartoon of the potential clump finding process in ID. 
For the potential and contours shown, the blue regions contain 
material which would be assigned to p-cores. 



of the fragment as it first becomes bound and the final mass 
of the sink particle that forms from it. Rather than breaking 
the cloud up into structures defined by density, we instead 
look at peaks in the gravitational potential. There are two 
advantages to such a method. First, the gravitational poten- 
tial distribution in the cloud is considerably smoother than 
the density distribution, since density fluctuations that do 
not carry sufficient mass cannot significantly contribute to 
potential field. Second, the strength of the gravitational po- 
tential determines whether a clump will collapse and how 
mass will flow. In a density distribution it is unclear which 
scale of the structure is important, whereas with potential 
the scale at which structures are bound is a clear physical 
quantity. Naturally, there is a disadvantage to this process: 
it becomes difficult to compare our structures to those ob- 
served in molecular clouds. However we believe some prop- 
erties may be comparable between the observed objects and 
those which are extracted from the simulation in this study. 
We discuss these in Section [U 

Potential wells are identified using an algorithm similar 

1994"), which works di- 



to CLUMPFIND (Wilhams et al 



rectly on SPH particles (|Klessen & Burkert 2000). We have 
modified the algorithm to use gravitational potential to iden- 
tify clumps rather than emission or density. In this scheme 
the SPH particle with the deepest gravitational potential 
forms the head of a clump, then the particle with the next 
deepest potential is either assigned to the same clump if it is 
a neighbour, or forms a new clump if not, and so on. Clumps 
are defined down to either a minimum positive potential, 
or the lowest contour which it shares with a neighbouring 
clump. Unlike the traditional CLUMPFIND algorithm, we 
use contour levels primarily to define the level at which po- 
tential clumps join, rather than to distinguish clumps from 
noise, and so our contours are numerous and finely spaced. 
This has the effect of subtracting the background potential. 
This is necessary as gravity is a long range force, and so is af- 
fected by both the mass inside and surrounding the p-core; 
hence we must remove the background to obtain the net 
effect on the mass within. P-cores, therefore, represent the 
local maximum above the surrounding background. Figure |2] 
shows a simple ID cartoon of the potential clump finding al- 
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gorithm where three potential clumps have been identified. 
Figure [S] shows an example of the potential clumps found 
from the central region of Figure [l] 

It is now wise to define our future terminology. Usually 
(though not always) a 'clump' is said to contain 50 — 500 
M© within 0.3 — 3 pc, and a 'core' is said to contain 0.5 — 5 
Mo within 0.03 - 0.2 pc ( [Bergin k Tafalla|2007t . Following 
this convention, as the structures identified by the potential 
clump find are typically small, we shall call them 'p-cores'. 
Note that there is no intrinsic requirement that a p-core 
is bound, as it may have sufficient internal energy (both 
thermal and kinetic) to prevent collapse. Once our p-cores 
form sinks, we only trace the total sink mass formed from 
them, because the spike in potential due to the sink particle 
distorts the boundaries; therefore all our p-cores are pre- 
stellar. 

The major aim of this paper is to investigate the link 
between pre-stellar p-cores and the stars formed from them. 
In order to do this we use two data-sets, the first of which is 
the 'composite' data-set. P-cores are found at 0.1 dynamical 
time {tdyn — 4.7 x 10^ yrs) intervals between 0.6 tdyn when 
star formation had just begun and 1.4 tdyn when the simula- 
tion was terminated. Snapshots from each time are combined 
into one data- set, which removes any time dependent effects 
and increases our data-set. The cores in the composite data- 
set have an average dynamical time of 7.4 x 10^ yrs meaning 
that a few long lived, non-transient cores will be included 
twice, but at different evolutionary states. In total there are 
573 p-cores. The composite data-set represents the view of 
cores in a molecular cloud at single points in time. It con- 
tains p-cores at different evolutionary states, with different 
levels of binding, many of which are transient. 

The second data-set is the 'bound' data- set, which con- 
tains the details of the p-cores at the point they first became 
bound. The p-cores are traced throughout the lifetime of 
the simulation, and if more than 80% of the mass belonging 
to a p-core remains grouped together in the next simula- 
tion time-step, the p-core is said to survive. The binding 
is traced throughout the p-core lifetime, and at the point 
where Erat > 1 for the first time its properties are recorded. 
Erat is defined as: 



(3) 



. + Ek 



where Ek is the kinetic energy calculated with respect to the 
center of velocity of the clump, Etherm is the thermal energy 
of the clump and Ep is the potential energy of the clump cal- 
culated using the relative depth of the potential well once 
the background has been subtracted. Hence the clumps iden- 
tified here are bound with respect to the environment in 
which they are formed, not merely when considered in isola- 
tion. This means that tidal forces from surrounding cores are 
taken into account when determining whether the p-core is 
bound. As the bound p-cores are identified throughout the 
simulation this data-set is also time independent. Due to 
this data-set being a synthesis over time and the p-cores be- 
ing identified using a quantity impossible to observe, these 
objects could not be found observationally. Despite this, the 
analysis is worthwhile as the bound data-set allows us to 
identify the fragmentation scale of the Molecular Cloud. The 
sink particles form when these bound cores collapse. In total 
there are 306 bound cores in our simulation. 
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Figure 4. The density of the sph particles assigned to a typical 
p-core plotted radially outward from the peak of gravitational 
potential. There is considerable dispersion due to substructure, 
but there is a clear trend showing a flattened central peak and 
density decreasing outwards. 



4 PHYSICAL PROPERTIES OF THE 
POTENTIAL CORES 

Unfortunately, a potential core can only be identified when 
positions and velocities are known in three dimensions, 
which is impossible observationally. Therefore, we now cal- 
culate the properties of the p-cores to allow a comparison 
to observational data. The average properties are shown in 
Table H 



4.1 P-core Shapes 

Figure [3] shows a close up view of the p-cores identified at 
a snapshot in time in one of the clustered regions of our 
molecular cloud. The p-cores are only quasi-spherical and 
often are elongated due to the filaments they are formed in. 
Within the p-cores there is still significant density substruc- 
ture, and further fragmentation is often observed after the 
p-core is first bound. Due to this substructure if we had used 
a traditional clumpfinding algorithm they would have been 
split into smaller cores. 

In order to quantify the extent of the p-cores, we use an 
effective radii, re//, within which 68% of the mass is con- 
tained. As the p-cores are not relaxed, we use the peak in the 
gravitational potential to define the centre rather than the 
centre of mass, which produces smoother density profiles. 
Figure |4] shows the density of the sph particles assigned to 
a typical p-core, plotted against their distance, r, from the 
central potential peak. The large dispersion shown in the 
density profile is due to the substructure within the core 
and the effect of the cloud being embedded in a non-uniform 
medium. 

Nonetheless despite the scatter there is a trend towards 
a fiattened density profile in the central region, as has been 
seen in observations (e.g Ward-Thompson et al.||T994 ). To 
get an idea of how centrally concentrated the p-cores were, 
a power law of the form p oc was fitted to the p-core and 
the best fit value of n determined. We excluded the radius 
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Figure 3. A close up view of some of the pre-stellar pcores identified from Figure^ The left panel shows a zoom of a clustered region in 
column density and the right panel shows the pre-stellar p-cores identified in it. Colours depict column density and the scale runs from 
0.001 gcm~^ (blue) to 10 gcm~^ (yellow). The location of the potential peak of the twenty p-cores in this region are shown by hollow 
circles. The p-cores are only quasi-spherical and exhibit significant substructure. 



Table 2. The average clump properties of the p-cores in the 
bound and composite datasets. Re// is the radius within which 
68% of the mass is contained. The density profile is the best fit 
value of n for the profile p (x. 





Bound 


Composite 


Mass (M©) 


0.70 


0.78 


Re// (AU) 


2.4 X 10^ 


3.7 X 10^ 




0.27 


0.40 


Dynamical Time (yrs) 


2.1 X 10^ 


7.4 X 10^ 


Density Profile (n) 


1.37 


1.08 



within which the first 10 % of the mass was contained from 
our fit due to the aforementioned flattening. 

Due to the dispersion from the density substructure and 
the core's non-spherical nature, there was a large degree of 
uncertainty, often as large as 50%. Despite this uncertainty, 
when plotted as a histogram Figure [s] (c), we find that the 
bound cores have exponents clustered around n = 1.36 ± 
0.35. This is intermediate between a shallow density increase 
(n = 1) and that of a free-falling envelope (n = 1.5), and is 
an indication that our p-cores are still being formed when 
they become bound. Typically exponents of n ^ 1.6 are 
expected for Class and Class 1 YSO's ( [Young et al.|20Q3 ) 
. Conversely, the profiles of the composite p-cores have a 
wide range of density exponents with an average value of 
1.08. 



4.2 Masses and Sizes 

Figure [s] (a) shows the clump mass function (CMF) of the 
bound and composite p-cores, in which there is a clear re- 
semblance to the IMF. Our previous work (Smith et al. 



ways obtained from MC structure using the traditional 
CLUMPFIND algorithm, but it was unclear whether this 
had any physical meaning. We demonstrate here that the 
bound cores also follow this distribution. This shows that at 
some level there is a link between molecular cloud structure, 
the formation of bound cores, and ultimately star formation. 
We examine the link between p-cores and stars in Section |5] 
The average p-core masses for the bound and compos- 
ite data-sets are very similar; 0.70 and 0.78 Mq's respec- 
tively. This is broadly consistent with the characteristic stel- 



2008t 



showed that an IMF-like mass distribution was al- 



lar mass (Chabrier 2003| . Both data-sets show this distri- 
bution, regardless the fact that members of the composite 
population are often unbound and do not form stars. 

The distribution of the radii of the p-cores is shown in 
Figure [5] The bound p-cores have effective radii in the re- 
gion of 2.4 X 10^ AU, and the composite p-cores have radii 
of about 3.7 x 10^ AU. Both distributions resemble a log- 
normal. Unlike the clump mass function, the distribution 
of the bound and composite data-sets differs in magnitude, 
which is due to the composite population containing a large 
number of unbound diffuse clumps. 



4.3 Binding 

The binding of the composite data-set is shown in Figure |6] 
(a), only 24% of the p-cores are bound {Erat > !)• Moreover, 
for these p-cores, there is a tail which includes highly bound 
objects. The composite population consists of p-cores in a 
wide range of evolutionary states, a large fraction of which 
will not go on to form stars. Figure |6] (b) also shows the 
p-core masses plotted against energy ratio. There is no cor- 
relation between binding and mass in our simulation, con- 
trary to some observations where the most massive cores 
are the most bound (e.g. Lada et al.||2008" ). This is due to 
the tendency of our larger cores to be more diffuse. Due 



6 Smithy Clark & Bonnell 




1 10 
Potential Well Mass [M^^,] 




1000 10" 
Effective Radius [AU] 




0.5 1 1.5 2 

Density Exponent [n] 



Figure 5. Histograms of the masses & sizes of the bound (solid line) and composite (dotted line) p-cores. Panel (a) shows the clump 
mass function. Masses above 0.2 Mq are resolved and the Saltpeter slope is denoted by a dashed line. Panel (b) shows the effective radii. 
Panel (c) shows the best fit values of n for the profile p oc r~^. The p-core mass function resembles the stellar IMF and the p-cores are 
typically small, centrally concentrated objects. 
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Figure 6. The binding of the composite p-cores. Panel (a) shows a histogram of the energy ratio of the cores, Erat ^ 1 are bound, where 
Erat = \Ep\/ Etherm + ^fc • Panel (b) shows the p-core masses plotted against energy ratio, blue circles denote cores with a steep n > 1.5 
density profile, green squares intermediate 1 < n < 1.5 profiles, and red crosses shallow profiles. Panel (c) shows the density exponent 
n plotted against energy ratio with error bars due to the poor fit from density substructure and non-spherical core shapes, the straight 
line has a gradient of two. There is no correlation between binding and mass, but there is a link to central concentration. 



to binding being roughly constant with p-core mass, bound 
cores appear to be selected with a uniform probability from 
the composite p-core distribution, which explains the resem- 
blance between the composite and bound CMF's. 

The density profile of the mass assigned to the core ra- 
dially outwards, on the other hand, does show a correlation 
with binding. In Figure [g] (b) , shallow (n < 1), intermedi- 
ate (1 < n < 1.5) and steep (n > 1.5) density profiles are 
denoted by circles, squares and crosses respectively. We find 
that the steepest profiles are nearly all bound, and about half 
of the intermediate profiles are bound. Figure |6] (c) shows 
the density exponents of the p-cores plotted against their 
energy ratio. When the profiles are shallow, a core is almost 
never bound, but when the potential well deepens and the 
core becomes more bound, the density profile steepens as the 
p-core becomes more centrally concentrated. However, there 
is a high degree of dispersion in this relation due to the diffi- 
culty of obtaining density fits from only quasi-spherical cores 



which contain substructure. Nonetheless, the steeper density 
profiles and increased central condensation show the clear- 
est indication of binding from all the observationally visi- 
ble quanities. This echoes the Bonnor-Ebert fBonnor||1956 



Ebert|1955) sphere model often used to model core observa- 



tions (e.g. Johnstone et al.|[2000p , in which collapse begins 



above a critical ratio of maximum central density to mean 
density. However our objects are not generally in equilib- 
rium, but are dynamically evolving in a similar manner to 
that shown in Ballesteros-Paredes et al. (2003). 



Figure [7| shows the internal ID velocity dispersions of 
our bound and composite p-cores. Originally 3D velocity dis- 
persions were calculated but these have been converted to 
the ID for easier comparison to observations. The mean 3D 
velocity dispersions were 0.27 and 0.4 kms~^ respectively, 
and when converted to a ID velocity dispersion this be- 
comes 0.16 and 0.23 kms~^. The sound speed of an isother- 
mal gas at lOK is 0.2 kms~^, meaning the potential cores 
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Figure 7. The one dimensional internal velocity dispersions of 
the p-cores in the hound(solid line) and composite (dotted line) 
datasets. The sound speed of an isothermal gas at lOK is 0.2 
kms~^ which means our p-cores are generally subsonic 



Figure 8. The cumulative mass functions from snapshots at solid 
0-6 tdyn, dotted 0.8 tdymshort dashed 1 tdyn and long dashed 1.2 
tdyn- The dot-dashed line shows the Salpeter slope. The mass 
function gets steeper with time as the high mass p-cores are 
formed earlier. 



are typically just on the verge of becoming supersonic, and 
are therefore coherent objects ( Goodman et al._1998 ). 

The mean dynamical timescale of a bound cores is cal- 
culated using the depth of the potential, rather than the 
standard practice of using the density, as the gravitational 
potential is the more clearly defined quantity for our p-cores. 
We equate the potential and kinetic energies using the virial 
theorem to obtain a typical velocity which is used to find the 
dynamical time as shown below, where Reff is the effective 
radius of the core and the gravitational potential. 



tdyn 



R. 



eff 



(4) 



The average dynamical time is 2.8 x 10^ yrs, and for the com- 
posite data-set it is about two times larger. This is more than 
ten times shorter than the dynamical time of the molecular 
cloud as a whole; meaning that several generations of cores 
can form and evolve throughout the lifetime of the simulated 
molecular cloud. 



4.4 The Core Mass Function with Time 

Figurejs] shows the mass functions of the pre-stellar (without 
sinks) p-core snapshots before they were merged into the 
composite data-set to integrate out time effects. Snapshots 
are shown at 0.6 to 1.2 simulation dynamical times (tdyn = 
4.7 X 10^ yrs) at intervals of 0.2 tdyn • The high mass end 
of the clump mass function gets steeper with time, meaning 
that more massive cores were formed at the beginning of the 
simulation than at the end. This is probably due to cores at 
later times being more likely to form in the vicinity of an 
existing potential well. This will tidally truncate the size of 
the core, and hence limit the mass it can accrete. In our 
simulation the major regions of star formation are formed 
at the same time, and hence their evolution is artificially 
synchronous. If this were not the case the steepening of the 
mass function would not have been visible. We aim to further 
discuss the effect of time on the mass function in a future 
paper. 



5 CLUMP MASSES & STELLAR MASSES 

Since the first observations of dense cores, a direct causal 
relation between them and stars was proposed (e.g. |Motte| 



et al. 1998). Alves et al. (2007) go further, and find that 



core masses and stellar masses are related with an efficiency 
of 0.3. However, this neglects the influence of environmental 
factors on the core during the accretion process. Moreover, 
Swift & Williamsl (|2008|) have shown that cores can have 



variable efficiencies or form multiple stars and still generate 
the expected IMF. 

We now examine the correlation between core masses 
and their resulting stellar masses. To investigate how much 
of the p-core mass is available to form stars, we trace the 
mass of the sink particles formed from them with time. If a 
p-core forms more than one sink, we add their masses. This 
means we are actually tracing the correspondence between 
p-cores and the stellar systems formed from them. Figure |9] 
shows the masses of the sink particles formed as a function 
of their p-core masses evaluated at a single snapshot in time. 
Sink masses are recorded at the end of the simulation and 
the p-cores in the snapshot population that did not form 
sinks are neglected. 

There is no clear relation between the p-core masses and 
the sink masses. The sink masses all lie near or well above 
the p-core mass values showing the importance of subse- 
quent accretion. Note that as feedback is not included and 
all sinks that form from a core are counted, it would be 
surprising if the efficiency of sink formation from cores was 
much less than 100%. Instead it is more than a factor of two 
higher. The actual stellar mass formed would of course be 
less than the sink mass, which simply represents the mass 
which would reach the inner disk of an accreting proto-star. 
During the accretion process mass will be ejected via a jet, 
a process which is not modelled here. However, if there is 
no correlation between the mass reaching the forming proto- 
star there is unlikely to be a correlation with the final star 
either. In fact the correlation is most likely to be further 
worsened by this process. For instance |Matzner McKee| 
( |2000. ) find a star formation efficiency of between 25% — 70% 
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Figure 9. The connection between p-core masses at a snapshot 
in time and their sink mass when the simulation was stopped. The 

sohd hne shows a 1-1 correspondence. There is a poor correlation Figure 11. The CMF of the p-cores with mass bins denoted by 

between p-core mass and the total sink mass formed from them. different colours. 



in model cores with feedback, depending on geometry. Inter- 
estingly, they also find that the stellar IMF formed from the 
cores is insensitive to core efficiency. Since all observations 
are necessarily at a snapshot in time this work suggests that 
the final stellar mass cannot be predicted accurately from 
the observed core mass even if our p-cores were directly ob- 
servable. 

One complication of Figure [5] is that the cores are all 
at different evolutionary states. We can address this by us- 
ing their properties at the point in time where each is first 
bound, as this is when fragmentation and collapse begins, 
and when the best correlation to the p-core mass would be 
expected. We record the total mass in sinks formed from 
these p-cores at successive core dynamical times after they 
are first bound. In our simulation there is no mechanism to 
halt accretion, but as the p-core dynamical times ( Equation 
|4| are short, we should be finding the sink masses before this 
is an issue. Figure shows the masses of the sinks after 1, 
2, 3 & 5 dynamical times. 

When comparing the sink mass at equal evolutionary 
stages of the bound p-cores, the correspondence is a lot 
tighter than that in Figure |9] At 1 dynamical time after 
the p-core is bound the sinks are below the 1-1 correspon- 
dence line. At 2 dynamical times they are around this level. 
These two stages, therefore, follow the accretion of the initial 
bound core, and so it would be surprising if this correlation 
did not exist! However, there is considerable scatter in this 
trend, with sinks masses being as much as a factor of two 
away from the trend in either direction. This indicates that 
some cores are accreting their envelope at a faster rate in 
these initial stages. 

As suspected, therefore, there must be a direct connec- 
tion between bound cores and star formation. This improved 
correlation of p-core mass and sink mass is found at a specific 
point in time, when first bound. This leads us back to the 
original scheme of star formation proposed by 'Jeans" ("1902), 
as our potential bound cores represent the local Jeans mass 
at the point of collapse. 

At 3 and 5 dynamical times the cores are accreting from 
their wider environment i.e. material which became bound 
to the cores subsequent to when Erat > 1 for the first time. 



The correlation between core masses and sink masses be- 
comes increasingly dispersed in log-space as accretion wipes 
out the direct correspondence and increases the offset as the 
sink-masses grow with time. Cores surrounded by a plentiful 
reservoir of gas on all sides will be more successful at these 
stages than those in a narrow filament for example. 

Despite the poor correlation between core mass and sink 
mass, a Salpeter like mass function is always maintained. 
We illustrate this graphically in Figure [TT] which shows the 
cumulative clump mass function of the bound cores. Each 
point on the graph represents a p-core mass when first bound 
(our best case scenario) and different mass ranges are de- 
noted with different colours. We now trace the evolution of 
these cores into sinks, keeping the colours of the sink the 
same as its parent core. If there was a perfect 1 — 1 correla- 
tion, then the colour bands would remain distinct. Figure [12] 
shows the resulting mass functions at subsequent dynamical 
times. The colours are now well mixed, showing that the p- 
cores are evolving with variable efficiencies. Nonetheless, the 
shape of the IMF is maintained throughout, due the the ef- 
fects of competitive accretion. Therefore, we could say that 
for a population of cores there is a high probability that a 
more massive core will form a more massive star, but for a 
specific object no reliable predictions of final mass can be 
made. 

In summary, the bound p-cores all form stars, but some 
are more successful at accreting additional mass than oth- 
ers. This is most likely due to environmental factors, for ex- 
ample, core geometry, surrounding gas reservoir, dynamical 
interactions or competition from neighbouring proto-stars. 
This suggests, that in order to understand star-formation, 
not only the cores of gas must be studied, but also the wider 
environment in which they form. The correspondence be- 
tween cores and stars is even more difficult to determine at 
a snapshot in time when cores are at a variety of evolution- 
ary stages. 
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Figure 10. The connection between clump mass and sink mass at successive dynamical times. Panel (a) tdyn = I5 (b) tdyn = 2, (c) 
tdyn = 3, (d) tdyn = 5. The solid line shows a 1-1 correspondence. There is now a clear connection between p-core mass and sink mass, 
but it still shows significant dispersion. 



6 DISCUSSION 



It has been theorised that the power law clump mass func- 



tion is formed as a product of supersonic turbulence ( |Hen- 

L,:i I h r\or:} l I t r>r>o l I tth I h nni-v l \ t^i "i or>r>i 1 



riksen|[T986l [Lar son" 1992'; 'Elmegreen' 1 19971 |Kle ssen"'2QQl' 
Padoan & Nordlund 2002; Hennebelle & Chabrier 2008 



The core mass function has been observed by many authors 
(e.g. [Motte et al.||1998| [Testi Sargent,, 1998^ [Johnstone 



et al. I [20001 [Nutter &; Ward-Thompson[|2007| jlkeda et al 
2007) and bears a remarkable resemblance to the stellar 



IMF. We now confirm that this profile is also observed in 
the first bound fragments we have identified in our simula- 
tion. 

However, the p-cores are not directly comparable to 
observed cores as they are identified in three dimensions, 
and use gravitational potential instead of density. Further, 
it has also been shown that the core properties obtained 
from applying CLUMPFIND to 2D observations are unreli- 
able (Kainu lainen et al.|2009] [Smith et al.|2008t . Our p-cores 
are more well defined as they are found from a smooth dis- 
tribution where all structures are significant. Comparisons 
between observed core properties and p-core properties must 
therefore be made cautiously. Nonetheless it is useful to con- 



trast the typical sizes of observed cores to identify which 
objects are on the most similar size scales to our first bound 
fragments. Table [S] summarises observed properties of den- 
sity cores in different regions and species from a variety of 
authors. 

P-cores are typically smaller than most density cores 
identified observationally, having average masses of 0.7 Mq 
and radii of 2.4 x 10^ AU when first bound. Their inter- 
nal ID velocity dispersions are typically just subsonic. The 
p-cores are most similar to high resolution observations of 
cores in nearby molecular clouds (e.g. Simpson et al. 2008| 
Enoch et al.[|2008 ). We expect that the difference between 



most observed cores and our potential cores is mainly due 
to resolution. Most surveys, particularly of more distant re- 
gions such as Orion, would not be able to resolve cores of 



the size found here, and as we showed in Smith et al. (2008) 



coarse resolution observations still produce a Salpeter like 
mass function, but with systematically higher masses. 

Although the clump mass function observed from den- 
sity fluctuations is shown to resemble the stellar initial mass 
function, there was no requirement for this to be true for the 
distribution of masses at which cores first become bound. In 
fact, in the competitive accretion theory of star formation 
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Figure 12. The initial mass functions of the sinks formed from bound cores the colours are the same as their parent p-cores shown in 
Figure [TTl The IMF's are recorded at intervals of 1,2,3 &; 5 p-core dynamical times. The mass bins get mixed up with time, but the 
shape of the mass function always resembles the stellar mass function despite the inexact correspondence between core mass and sink 
mass. 



Table 3. A comparison to a small selection of core measurements from the literature. Authors are shown, additional velocities are from 
* |Lada et al .' (2008), *'^Rosolow sky et al.| ( |2008| and show the ID internal velocity dispersion. [Simpson et al~| (2008 ) is a re-analyses of 
mm observations in Ophiuchus from various published works. Note the p-cores found in this work are defined differently, as they use 3D 
potential rather than 2D column density. 



Author 



Species 



Typical Clump Size Clump Mass Range Mean Velocity Dispersion 







[AU] 


Mo 


kms"*^ 


OrionA 


H13CO+ 


2.9 X 10^ 


2.1-81 


0.52 


Orion B 


mm cont 


1.4 X 10^ 


0.4 - 30 




Ophiuchus 


extinction 


- 5 X 10^ 


0.5 - 28 


< 0.2* 


Ophiuchus 


mm 850mic 


~ 3 X 10^ 


0.1 - 10 




Oph/Per/Serp 


1.1mm cont 


- 7 X 103 


0.2-4.8 


0.24** 


simulated 


potential (bound) 


2.4 X 10^ 


0.2 - 5 


0.16 


simulated 


potential (composite) 


3.7 X 10^ 


0.2 - 10 


0.23 



Ikeda et al. 


(20071) 


Johnstone e 


t al. 


j2006fc 


Xlves et al. 


(20071 


STmpson et al. (2008) 


^Tnoch et al. 


,(,2008^ 



This work 
This work 



the seeds of gravitational collapse can follow any distribu- 
tion and the IMF will still be generated through subsequent 
accretion ( Zinnecker|1982 |. That the time integrated bound 
fragments also follow this distribution further supports the 
idea that they are sampled from a hierarchical density dis- 



tribution generated by turbulence. There is some indication 
that this may not be true at late snapshots in time. Few large 
diffuse p-cores are formed later in the simulation, which re- 
sults in a steeper mass function at late epochs. 

However, at a snapshot in time the majority (76%) of 
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the p-cores are unbound and may be transitory. This has 
imphcations for observational surveys which are necessarily 
at a snapshot in time, as the vast majority of objects have 
either not reached the point of collapse or will not form stars 
at all. The case for a 1-1 correlation between core mass and 
stellar mass can not hold in these cases. There is also the fur- 
ther complication that our p-cores contain substructure that 
could be further identified as individual objects if one were 
performing a density decomposition (depending on the res- 
olution of the data). Naturally, these yet smaller objects are 
even less bound than their parent p-core. The high value of 
unbound pre-stellar p-cores in our analysis is surprising be- 
cause as we identify them using gravitational potential this 
force must be significant. There is observational evidence for 



unbound cores. Lada et al. (2008) find that the majority of 
cores found from dust extinction in the Pipe nebula are un- 
bound, and Andre et al. (2007) find 25% of their cores in 



Ophiuchus to be unbound. 

The p-cores showed significant density sub- structure, 
despite being identified from a smooth potential well. If a 
traditional clump finding algorithm using emission or den- 
sity had been used to identify them, they would have been 
split into smaller objects. The p-cores were also only quasi- 
spherical, and their central peak in potential did not neces- 
sarily correspond with their center of mass, as they were not 
relaxed. This suggests that using smooth symmetrical mod- 
els to study core collapse (such as Bonnor-Ebert spheres) 
can at best only be a first approximation. 

Unfortunately, our simulation does not include mag- 
netic fields, which can act against collapse on some length 
scales ( jHennebelle Passot||2006| ). [Price Bate] (|2008f 
found that magnetic fields introduce voids into molecular 
cloud structure and reduce subsequent accretion. We also 
do not have radiative transfer in our simulation. However, 
as we are considering the cores early in their evolution before 
a protostar is formed, and accretion luminosity only affects 
scales of a few 100 AU ( [Krumholz et al.||2007| |Bate||2009| ), 
this should not be a significant issue. 

Due to the similarities between the Stellar IMF and the 
clump mass function it has been proposed that the mass of 
stellar systems is directly related to the clump they formed 
from (e.g. [Motte et al.||1998] |Alves et al.||2007t . However, 
when we traced the total mass in sink particles formed from 
our p-cores at a snapshot in time, as would be visible obser- 
vationally, the correspondence was effectively non-existent. 
A stronger correlation between core masses and total sink 
mass was found when the p-core masses were recorded at the 
time when they were first bound. Further accretion increases 
the dispersion in this relation and introduces an offset. 

Throughout the simulation the most massive bound 
core had a mass of only 6.35 M©, whereas the most mas- 
sive sink at the end of the simulation had a mass of 27.97 
M0. We did not find any p-core which could form a massive 
star simply from its own material when first bound. All the 
bound cores were at best only marginally supersonic and 
hence were not supported by turbulence, contrary to the 



massive turbulent cores predicted by McKee & Tan ( 2003 ) 



Note, however, that the lack of feedback and magnetic fields 
in these simulations may limit the ability to sustain turbu- 
lence in cores. 

In most cases gas from outside the region initially bound 
was accreted by the sinks formed from the p-core after only 



6 X 10^ yrs. This subsequent accretion means that the en- 
vironment surrounding the p-core is also important for its 
future evolution, for example whether it is surrounded by a 
large gas reservoir or is in a narrow filament. 

There is a higher probability that a more massive p-core 
will form a larger total sink mass, but for an individual p- 
core no accurate prediction can be made. The link between 
core mass and sink mass is poor because there are a num- 
ber of additional environmental factors beyond the bound 
mass which will affect the evolution of the core, for example 
the shape of the core and the distribution of gas it is embed- 
ded within. Despite the poor correspondence between p-core 
masses and sinks, the shape of the mass function obtained 
from the p-cores as they evolve into sinks always resembles 
the stellar IMF. 



7 CONCLUSIONS 

We identified the earliest fragmentation in molecular clouds 
using the distribution of gravitational potential rather than 
density. This is a smoother distribution, due to gravity being 
a long range force, and it allows us to assign boundaries with 
clearer physical meanings. We called our fragments 'p-cores', 
and traced their evolution with time. If they became bound, 
we identified their mass when (Ep > Ek -\- Etherm) for the 
first time, and we related this to the mass in sink particles 
formed from them. This allowed us to identify the scale of 
the initial fragmentation in molecular clouds, and trace its 
evolution towards the stellar IMF. Our main conclusions are 
as follows: 

(i) The time-integrated mass function of just-bound grav- 
itational potential 'p-cores' resembles the stellar initial mass 
function in a similar manner to that of gas density cores 

(ii) The bound p-cores are most similar to the very small- 
est density cores currently observable ( an average of 0.7 M© 
with a radius of 2300 AU). 

(iii) P-cores exhibit significant density substructure. 

(iv) There is a poor correlation between p-core mass and 
the total mass in sinks formed from them. This is particu- 
larly true when the p-cores are recorded at a snap-shot in 
time, but still holds when the mass is recorded when they 
are first bound. 

(v) The sink particles formed from the p-cores accrete 
from beyond the region initially bound. This means that 
the surrounding environment of the core also has an effect 
on the star formation in our simulation. 

(vi) Despite the poor correlation between p-core mass and 
sink masses the sink IMF always resembles the stellar IMF. 
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